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We show how it is possible to suppress decoherence using tailored external forcing acting as 
pulses. In the limit of infinitely freqüent pulses decoherence and dissipation are completely frozen; 
however, a significant decoherence suppression is already obtained when the frequency of the pulses 
is of the order of the reservoir typical frequency scale. This method could be useful in particular to 
suppress the decoherence of the center-of-mass motion in ion traps. 
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' I. INTRODUCTION 

Decoherence is the process which limits our ability to maintain puré quantum states, or their linear superpositions. 
It is the phenomenon by which the classical world appears from the quantum one In more physical terms it is 
described as the ràpid destruction of the phase relation between two, or more, quantum states of a system caused 
by the entanglement of these states with different states of the environment. The present widespread interest in 
' decoherence is due to the fact that it is the main limiting factor for quantum information processing. We can store 
^ , information, indeed, in two-level quantum systems, known as qubits, which can become entangled each other, but 
■ decoherence can destroy any quantum superposition, reducing the system to a mixture of states, and the stored 
information is lost. For this reason decoherence control is now becoming a rapidly expanding field of investigation. 

In a series of previous papers ||-|5| we have faced the control of decoherence by actively modifying the system's 
dynamics through a feedback loop. This procedure turned out to be very effective, in principle 0, to slow down the 
decoherence of the only one experiment up to present, in which the decoherence of a mesoscopic superposition was 
Q\ ' detected. The main limiting aspect of this procedure is conncctcd with the need of a measurement. In order to do 
the feedback in the appropriate way, one has first to perform a measurement and then the result of this measurement 
can be used to operate the feedback. However, any physical measurement is subject to the limitation associated 
with a non-unit detection efficiency. We have shown || that with detection efficiency approaching unity the quantum 
superposition of states stored in a cavity can be protected against decoherence for many decoherence times tdec, where 
tdec is defincd as the cavity relaxation time divided by the average photon number j7j . 

We wish now to face the problem of eliminating the measurement in controlling the decoherence. We show, here, 
how it is possible to inhibit decoherence through the application of suitable open-loop control techniques to the system 
^ , of interest, that is, by using appropriately shaped time-varying control fields. To be more specific, decoherence can 
be inhibited by subjecting a system to a sequence of very freqüent parity kicks, i.e., pulses designed in such a way 
that their effect is equivalent to the application of the parity operator on the system. 

The paper is organized as follows: In section II the parity kicks method is presented in its generality and it is 
shown how decoherence and dissipation are completely frozen in the limit of infinitely freqüent pulses. In Section 
III the method is applied to the case of a damped harmònic oscillator, as for example, a given normal mode of a 
system of trapped ions. In Section IV the numerical results corresponding to this case are presented showing that 
a considerable decoherence suppression is obtained when the parity kicks repetition rate becomes comparable to the 
typical timescale of the environment. In Section V the possibility of applying this scheme to harness the decoherence 
of the center-of-mass motion in ion traps are discussed. 
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II. THE GENERAL IDEA 



Let us consider a genèric open system, described by the Hamiltonian 

H = H A + H B +H INT , (1) 

where Ha is the bare system Hamiltonian, Hb denotes the reservoir Hamiltonian and ÍÏ/tvt the Hamiltonian describing 
the interaction between the system of interest and the reservoir, which is responsible for dissipation and decoherence. 
We shall now show that if the Hamiltonian (Eh possesses appropriate symmetery properties with respect to parity, it 



1 



is possible to actively control dissipation (and the ensuing decoherence) by adding suitably tailored time-dependent 
external forcing acting on system variables only. The new Hamiltonian becomes 

Htot = H + H kick (t) , (2) 

where 

oo 

H Hck (t) = H k 0(t — T — n(T + r )) 6 ((n + 1)(T + r ) - í) , (3) 

n=0 

(0(t) is the usual step function) is a time-dependent, periòdic, system operator with period T + tq describing a train 
of pulses of duration tq separated by the time interval T. The stroboscopic time evolution is therefore given by the 
following evolution operator 

U(NT + Nt ) = [U kick (T )U (T)} N (4) 

where U k i C k{To) — exp{— ítq(H +H k )/h} describes the evolution during the external pulse and Uq(T) = exp{—iTH/h} 
gives the Standard evolution between the pulses. We now assume that the external pulse is so strong that it possible 
to neglect the Standard evolution during the pulse, H k 3> H, and then we assume that the pulse Hamiltonian H k and 
the pulse width tq can be chosen so as to satisfy the following parity kick condition 

U ktck (T )~e-K H » To =P, (5) 

where P is the system parity operator. 

It is now possible to sec that such a timc-dcpcndcnt modification of the system dynamics is able to perfectly protect 
the system dynamics and completely inhibit decoherence whenever the following general conditions are satisficd: 

PH A P = H A (6) 
PH INT P — —H INT , (7) 

that is, the system Hamiltonian is parity invariant and the interaction with the external environment anticommutes 
with the system parity operator. To be more specific, one has that in the ideal limit of continuous parity kicks, that 
is 

T + t -> 

N -► oo (8) 
t = N(T + r ) = const. 

the pulsed perturbation is able to eliminate completely the interaction with the environment and therefore all the 
physical phenomena associated with it, i.e., energy dissipation, diffusion and decoherence. To see this it is sufhcient to 
consider the evolution operator during two successive parity kicks, which, using the parity kick condition of Eq. (||), 
can be written as 

U(2T + 2t ) = Pe-* HT Pe~^ HT ; (9) 

since 

Pe-* HT P = e -* PHPT = e -Í(HA+H B -H INT )T ^ ( 1Q ) 

one has that the time evolution after two successive kicks is driven by the unitary operator 

U(2T + 2t ) = e -ïï( H A+H B -H INT )T e --ï : (H A +H B +H 1NT )T _ ^ 

This expression clearly shows how the pulsed perturbation is able to "freeze" the dissipative interaction with the 
environment: the application of two successive parity kicks alternatively changes the sign of the interaction Hamilto- 
nian between system and reservoir. Therefore one expects that, in the limit of continuous kicks, this sign inversion 
becomes infinitely fast and the interaction with the environment averages exactly to zero. 

This fact can be easily shown just using the defmition of continuous kicks limit. In fact, in this limit 



U(t) = lim e -i(HA+H B -H INT )T e -±(H A +H B +H INT )T 2r + 2 -o 
T+TQ-+0 L J 
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which, using just the definition of the exponential operator, yields 

U(t) = e -UH A +H B )t 

This means that in the ideal limit of continuous parity kicks, the interaction with the environment is completely 
eliminated and only the free uncoupled evolution is left. 

Let us briefiy discuss the physical interpretation of this result. The continuous kicks limit (0) is formally analogous 
to the continuous measurement limit usually considered in the quantum Zeno effect (see for example HH) in which a 
stimulated two-level transition is inhibited by a sufficiently freqüent sequence of làser pulses. However this is only a 
mathematical analogy because at the physical level one has two opposite situations. In fact, during the pulses, in the 
quantum Zeno effect the interaction with the environment (i.e., the measurement apparatus) prevails over the infernal 
dynamics, while in the present situation the interaction with the reservoir is practically turned off by the externally 
controlled infernal dynamics (see Eq. (|^)). 

This parity kick idea has instead a strong relationship with spin echoes phenomena |Ï0| in which the application of 
appropriate rf-pulses is able to eliminate much of the dephasing in nuclear magnètic resonance spectroscopy experi- 
ments. These clever rf-pulses realize a sort of time-reversal and a similar situation takes place here in the case of parity 
kicks. In fact the relevant evolution operator is given by the unitary operator governing the time evolution after two 
pulses U(2T + 2ro) of Eq. (|IÏ|), in which the Standard evolution during a time interval T is followed by an evolution 
in which the interaction with the environment is time-reversed for the time T. If the time interval T is comparable 
to the timescale at which dissipative phenomena take place, this time-reversal is too late to yield appreciable effects, 
while if T is sufficiently small the very freqüent reversal of the interaction with the reservoir may give an effective 
freezing of any dissipative phenomena. 

Two very recent papers Jïl"[ have considered a generalization of this parity kicks method and have showed that 
a complete decoupling between system and environment can be obtained for a genèric system if one considers an 
appropriate sequence of infinitcly freqüent kicks, realizing a symmetrization of the evolution with respect to a given 
group. The present parity kicks method is in fact equivalent to symmetrize the time evolution of the open system 
with respect to the group , composed by the identity and the system parity operator P. Here we focus only on this 
case because the experimental realization of the parity kicks of Eq. (^) is quite easy in many cases and moreover the 
applicability conditions of the present method, Eqs. (g) and (Q), are satisfied by many interesting physical systems. 

The continuous kicks limit is only a mathematical idealization of no practical interest. However it shows that there 
is in principle no limitation in the decoherence suppression one can achieve using parity kicks. Moreover it clcarly 
shows that one has to use the most possible freqüent pulses in order to achieve a significative inhibition of decoherence. 
Therefore the relevant question from a physical point of view is to determino at which vàlues of the period T + tq one 
begins to have a significant suppression of dissipation and decoherence. We shall answer this question by considering 
in the next section, a specific example of experimental interest, that is, a damped harmònic oscillator, representing 
for instance a normal mode of a system of trapped ions, or an electromagnètic mode in a cavity. 



III. PARITY KICKS FOR A DAMPED HARMÒNIC OSCILLATOR 



Let us consider a harmònic oscillator with bare Hamiltonian 

H A = hu} a í a , (14) 

describing for example a given normal mode of a system of ions in a Paul trap or an electromagnètic mode in a 
cavity. In these two cases the relevant environmental degrees of freedom can be described in terms of a collection of 
independent bosonic modes p2[ 

H B =J2 t^kb{b k , (15) 

k 

representing the elementary excitations of the environment (in the cavity mode case they are simply the vacuum 
electromagnètic modes). Moreover, the interaction with the environment is usually well described by the following 
bilinear term in which the "counter-rotating" terms are dropped HQ 

H INT = Y,hlk{abl + a)b k ^ . (16) 

k 

In these cases time evolution is usually described in the frame rotating at the bare oscillation frequeney ujq in which 
the effective total Hamiltonian of Eq. <m) becomes 
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H = H' B + H INT (17) 

where 

H' B = Y J H"k - uJo)b{b k . (18) 

fe 

It is immediate to see that in this example the conditions of Eqs. @ and (0) for the application of the parity kicks 
method are satisfied. More generally, conditions (^) and (R) are satisfied whenever Ha is an even function of a and 
and í/jjvt is an odd function of a and . In the harmònic oscillator case the pulsed perturbation realizing the 
sequence of parity kicks can be simply obtained as a train of 7r-phase shifts, that is 

H klck {t)=E(t)a^a, (19) 

where E(t) = E Q J2n=o d {t — T — n(T + r )) 9 ((n + 1)(T + r ) — t) and the pulse height E Q and the pulse width r 
satisfy the condition 

Eqtq = (2n + l)7r?ï n integer . (20) 

To see how large the pulses repetition rate has to be in order to achieve a relevant decoherence suppression, first of 
all one has to compare the pulsing period T + tq with the relevant timescales of the problem. The harmònic oscillator 
dynamics is characterized by the free oscillation frequency lüq and by its energy decay rate 7 (which is a function 
of the couplings 7fc); the bath is instead characterized by its ultraviolet frequency cutoff w c which essentially fixes 
the response time of the reservoir and generally depends on the system and bath considered. A reservoir is usually 
much faster than the system of interest and this means that one usually has u c 3> 7. Typically the reduced dynamics 
of the system of interest is described in terms of effective master equations which are derived using the Markovian 
approximation (see for example ]Ï3| ) which means assuming the limit uj c — > 00; however, the existence of a finite cutoff 
lü c is always demanded on physical grounds [ fÏ5| and this parameter corresponds for example to the Debye frequency 
in the case of a phonon bath. 

From the preceding section, it is clear that the pulsing period T + tq has to be much smaller than I/7, otherwise the 
change of sign of the interaction Hamiltonian is realized when a significant transfer of energy and quantum coherence 
from the system into the environment has already taken place. At the same time it is easy to understand that the 
condition uj c (T + tq) <C 1 is a sufhcient condition for a significant suppression of dissipation, since in this case the sign 
of the interaction Himt changes with a rate much faster than every frequency of the bath oscillators; each bath degree 
of freedom becomes essentially decoupled from the system oscillator and there is no significant energy exchange, i.e. 
no dissipation. Therefore the relevant question is: for which value of the pulsing period T + tq within the range 
[I/7, l/uj c ] one begins to have a relevant inhibition of dissipation? 

We shall answer this question by studying in particular the time evolution of a Schròdinger cat state, that is, a 
linear superposition of two coherent states of the oscillator of interest 

\^ v )=N v (\a )+e^\-a Q )) , (21) 

where 

N v - 1 



^2 + 2e- 2 \ a o\ 2 costp 

We consider this particular example because these states are the paradigmàtic quantum states in which the progressive 
effects of decoherence and dissipation caused by the environment are well distinct and clearly visible ||. From the 
general proof of the preceding section it is evident that whenever one finds a significant suppression of decoherence 
in the Schròdinger cat case, this implies that the system-reservoir interaction is essentially averaged to zero and that 
one gets a significant system-environment decoupling in general. 

Describing the evolution of such a superposition state in the presence of the dissipative interaction with a reservoir 
of oscillators which is initially at thermal equilibrium at T — is quite simple. In fact, it is possible to use the fact that 
a tensor product of coherent states retains its form at all times when the evolution is generated by the Hamiltonian 
@, that is 

|«o) ® I] l^(°)> -> l a W) ® II !&(*)> . ( 22 ) 

fe k 

where the timc-dcpcndcnt coherent amplitudes satisfy the following set of linear differential equations 
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à(t) 



Pk(t) = -i(üJk - Lü )f3k{t) - ij k a(t) . 
Therefore, using Eq. (|22|), one has the following time evolution for the Schròdinger cat state ( pïl"| ) 

i^> ® n - ^ ( i«(*)> ® n i&(*)> + eiv i - «(*)> ® n i - & w> ) ■ 

fe V fe kJ 

The corresponding reduced density matrix of the oscillator of interest at time í, p(í), is given by 

p(t)=N^{\a(t)){a(t)\ + \-a(t))(-a(t)\ 

+D(t) [e-^\a(t))(-a(t)\+e^\-a(t))(a(t)\]} , 



(23) 



(24) 



(25) 



i.e. is completely determined by the coherent amplitude a(t) which is a decaying function describing dissipation of the 
oscillator 's energy and by the function D(t) describing the suppression of quantum interference terms, whose explicit 
expression is 



D(t) = Y[(Pk(t)\ - /3fc(í)> = exp l \Pk(t)\' 



(26) 



These results for the time evolution of the Schròdinger cat are vàlid both in the presence and in the absence of 
the external pulsed driving. In fact the only difference lies in the fact that in the absence of kicks one has the 
Standard evolution driven by exp {— iHt/h}, where H is given by (^7|), while in the presence of parity kicks one has 
a stroboscopic-like evolution driven by a unitary operator analogous to that of Eq. ( |ll| ) 



U{2NT + 2Nt ) 



'K H B T Op--k(H' B ~H RW A)T—-% : H B To p -j;(H B +H RWA )T 



N 



(27) 



(we have neglected only the system-reservoir interaction during the parity kick and therefore we have added two terms 
e -i H BTo m w i tn respect to Eq. Q). 

Due to the general result (|22|), it is convenient to express the state of the whole system in terms of a vector 
(a(£), • ■ ■ , (3k(t), ■ ■ ■) whose zero-th component is given by the amplitude a(t) and whose k-th component is given by 

the amplitude 0k(t), k = 1, 2, In this way the formal solution of the set of linear equations ( p3| ) can be expressed 

as 



(28) 



where ^({7^}, t) is a matrix whose first matrix element A({"/ k }, t)oo is given by the following inverse Laplace transform 

1 



/ "(0 \ 




(a \ 


Pk(t) 


= A({ lk },t) 





V ; J 




V i / 



where 



A({ lk },t) m = £- 



z + K(z) 



(29) 



z + i(uj k - ÜJq) 



(30) 



All the other matrix elements can be expressed in terms of this matrix element A({7fc}, í)oo in the following way: 



A({ lk },t) ok - A{{ lk },t) k{) = -i lk / dse"^— °)M({ 7fe },í- s) 



Hm 



(31) 



^({7fe},í)fefe' =6kk>e- ilu "·- Uo)t -Tk'Yk / dse-^*-»^-^ / ds' e^'-^-^ A{{ lk \, s ')c 

Jo Jo 



(32) 
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In the presence of kicks, the stroboscopic evolution is determined by the consecutive application of the evolution 
operator for the elementary cycle lasting the time interval 2T + 2tq (see Eq. (p7|)), during which one has a Standard 
evolution for a time T, followed by an uncoupled evolution for a time To, then a "reversed" evolution for a time T in 
which the interaction with the environment has the opposite sign and another uncoupled evolution for the time tq at 
the end. It is easy to understand that in terms of the vector of coherent amplitudes, the evolution in the presence of 
parity kicks can be described as 



/ a(2NT + 2Nt ) \ 
p k (2NT + 2N To 



\ 



[F(ro)A({- 7fc },T)F(ro)A({ 7fe },T)] 



í a °\ 


V ■ J 



(33) 



where F(tq) is the diagonal matrix associated to the free evolution driven by H' B and is given by F(tq)íj = 
Sí^j exp{— i(uji — u} )t }. We have seen that the quantities of interest are (see Eq. (|25|)) a(t) and D(t); for the 
amplitude of the two coherent states one has 



*(*) = A ({lk},t) 00 a 



in the absence of kicks and 



t(t) = {[F(ro)A({- lk },T)F(r )A({ lk },T)} N } oo a 



(34) 



(35) 



in the presence of kicks (t = 2NT + 2Ntq). As for the function D(t), it is clear that each amplitude /3fc(í) in ( |26| ) is 
proportional to ao and therefore one can write 



D(t) =cxp{-2|a | 2 r/(í)} 



where 



/(í) = El A (Í74,í) to 



feO 



(36) 
(37) 



in the absence of kicks and 



rt*) = E I { ({-Tfc}, T) F(r )A ({ 7fc }, T)] N \ I 2 



fcO 



(38) 



in the presence of kicks (t = 2NT + 2Nt ). 



IV. NUMERICAL RESULTS 



In the description of dissipative phenomena one always considers a continuum distribution of oscillator freqüències 
in order to obtain an irreversible transfer of energy from the system of interest into the reservoir. Moreover, most 
often, also the Markovian assumption is made which means assuming an infinitely fast bath with an infinite frequeney 
cutoff lü c . This case of a Standard vacuum bath in the Markovian limit is characterized by an infinite, continuous and 
flat distribution of couplings [Ofl 



which implies 



Using Eqs. @, © and @, we get 



7H 2 = 



K(z) = - 
V ' 2 



Q-Mark{t) = a 



T)Mc 



rfeW=E 



(39) 

(40) 

(41) 
(42) 



G 



The last sum in the expression for rjMark(t) has to be evaluated in the continuum limit, that is, replacing the sum 
with an integral over the whole real w-axis and using 7^ — * j/2tt. The result is the Standard vacuum bath expression 
for the decoherence function r)(t) [|| 

r)Mark(t) = 1 - e-T* . (43) 

From the general expressions of the above section it is clear that it is not possible to solve the dynamics in the 
presence of parity kicks in simple analytical form. We are therefore forced to solve numerically the problem, by 
simulating the continuous distribution of bath oscillators with a large but finite number of oscillators with closely 
spaced freqüències. To be more specific, we have considered a bath of 201 oscillators, with equally spaced freqüències, 
symmetrically distributed around the resonance frequency c^o, i.e. 

u k = ujo + fcA A = ^ (44) 

k max = ^ = 100 => Lü" k iax = 2íü (45) 

kmin = -L, = -100 =► ujr n = , (46) 

(47) 

and we have considered a constant distribution of couplings similar to that associated with the Markovian limit 

7Ï = ^ Vfc . (48) 

Considering a finite number of bath oscillators has two main effects with respect to the Standard case of a continuous 
Markovian bath. First of all, the adoption of a discrete frequency distribution with a fixed spacing A implies that 
all the dynamical quantities are periòdic with period T rev = 2n/ A flil which can be therefore considered a sort of 
"revival time" . It is clear that our numerical solution will correctly describe dissipative phenomena provided that 
we focus only on not too large times (say t < n/A). Moreover, the introduction of a finite cutoff (uj c = 2lüq in our 
case) implies a modification of the coupling spectrum "/(uj) at very high frequency with respect to the infinitely flat 
distribution of the Markovian treatment (see Eq. (|3^)). This fact will manifest itself in a slight modification of the 
exponential behavior shown by Eqs. (|4Ï]) and ( ^3| ) at very short times (í ~ w^ 1 ) |ïq ]. We have verified both facts in 
our simulations. However, to facilitate the comparison between the dynamics in the presence of parity kicks with the 
Standard case of a Markovian bath, we have choscn the paramcters of our simulation so that, within the time interval 
of interest, O.I/7 < t < 3/7 say, we found no appreciable difference between the Standard Markovian bath expressions 
( [ïl·| ) and (|4^) and the corresponding general expression for a discrete distribution of oscillators ([34|) and (|37j). 

Let us now see what is the effect of parity kicks on decoherence, by studying first of all the behavior of the decoherence 
function rj(t) for different vàlues of the pulsing period T + tq. As discussed above, we expect an increasing decoherence 
suppression for decreasing vàlues of T+tq and this is actually confirmed by Fig. 1, showing rj{t) for different vàlues of T 
(we have chosen tq — for simplicity in the numerical simulation). In fact the increase of rj(t) is monotonically slowed 
down as the time interval between two successive kicks T becomes smaller and smaller. Moreover Fig. 1 seems to 
suggest that decoherence inhibition becomes very efficient (r](t) ~ 0) when the kick frequency 1/T becomes comparable 
to the cutoff frequency uj c . To better clarify this important point we have plotted in Fig. 2 the value of the decoherence 
function at a fixed time as a function of the time between two successive kicks T. In Fig. 2a the decoherence function 
at half relaxation time i](t = O.5/7) is plotted as a function of ui c T/2ir, while in Fig. 2b the decoherence function 
at one relaxation time rj(t — 1/7) is plotted again as a function of u> c T/2tt. In both cases one sees a quite sharp 
transition at lu c T/2tt = 1: decoherence is almost completely inhibited as soon as the kick frequency 1/T becomes larger 
than the cutoff frequency u> c /2tt. Therefore we can conclude that by forcing the dynamics of an harmònic oscillator 
with appropriate parity kicks one is able to inhibit decoherence almost completely. This decoherence suppression 
becomes significant when the kick frequency 1/T becomes of the order of the typical reservoir timescale, i.e. the cutoff 
frequency lü c /2tt. 

A similar conclusion has been reached by Viola and Lloyd ]Ï7| , who have considered a model for decoherence 
suppression in a spin-boson model very similar in spirit to that presented here. In their paper they have considered 
a single 1/2-spin system coupled to an environment and they have shown how decoherence can be suppressed by a 
sequence of appropriately shaped pulses. They have shown that the 1/2-spin decoherence can be almost completely 
suppressed provided that the pulses repetition rate is at least comparable with the environment frequency cutoff. 
Similar results have been obtained in a more recent paper p8[ in which a generalization to more qubits and more 
general sequences of pulses is considered. However, despite the similarity of our approach to decoherence control to 
that of these papers, there are important differences between the present paper and Refs. ]Ï7j , |Ï8| |. First of all, we 



7 



have considered a harmònic oscillator instead of a 1/2-spin system, but above all, we have considered a dissipative 
bath, that is a zero temperature reservoir inducing not only decoherence (i.e. quantum information decay), but also 
dissipation (i.e. energy decay). In Ref. JÏ7t on the contrary, there is only decoherence and the 1/2-spin system 
energy is conserved. The more general nature of the present model hclps in clarifying a main point of this impulsivc 
mcthod to combat decoherence: parity kicks do not simply suppress decoherence but tend to completely inhibit any 
interaction between system and environment. This means that also energy dissipation is suppressed; in more intuitive 
terms, the external pulsed driving is perfectly "phase matched" to the system dynamics, so that any transfer of energy 
and quantum coherence between system and bath is inhibited. The fact that decoherence suppression in our model 
is always associated with the suppression of dissipation can be easily shown using the fact that the evolution of the 
whole system is unitary. In particular this implies that the matrix ^({7^}, í), driving the evolution in absence of kicks 
(see Eq. (28)), and the matrix [F(tq)A ({— 7fc}, T) F(tq)A ({7*,}, T)} , driving the evolution in the presence of kicks 
(see Eq. (33)), are unitary matrices and therefore subject to the condition 

|L o| 2 +^|L fco | 2 = l, (49) 



where Lij denotes any of the two above mentioned matrices. Now, using the reduced density matrix of Eq. fl25|), it is 
easy to derive the mean oscillator energy 

1 — ens í^p _2 l a °l 2 

(H A (t)) = ^ |«(í)| 2 — • (50) 

1 + cos (pe z \ a °\ 

From Eqs. ( |34| ) and ( [j5|) one therefore derives that the normalized oscillator mean energy is given by 

(Ha®) 



(Ha(0)) 



= l·4({7*},*)oor (51) 



{[F(r )A ({-7.}, T) F(r )A ({ 7fc }, T)f}J (52) 



in the absence of kicks and 

(H A {t)) 
(Ha(0)) 

in the presence of kicks (í = 2NT + 2Ntq). Now, by considering Eqs. ( |37| ) and ( |38| ) and the unitarity condition (|4 
it is immediate to get the following simple relation between decoherence and dissipation 

which is vàlid both in the presence and in the absence of kicks. This equation simply shows that when decoherence 
is suppressed (j](t) ~ 0) the oscillator energy is conserved. 

A qualitative demonstration of the ability of the kick method to suppress dissipation and decoherence is provided 
by Fig. 3. (a) shows the Wigner function of an initial odd cat state with cto = y5, tp = 7r; (b) shows the Wigner 
function of the same cat state evolved for a time t = 1/7 in the presence of parity kicks with u c T — 3.125 and (c) the 
Wigner function of the same state again after a time t — I/7, but evolved in absence of kicks. This elapsed time is 
ten times the decoherence time of the Schròdinger cat state, td ec = C^lWo] 2 )" 1 0, i.e, the lifetime of the interference 
terms in the cat state density matrix in the presence of the usual vacuum damping. As it is shown by (c), the cat 
state has begun to loose its energy and has completely lost the oscillating part of the Wigner function associated to 
quantum interference. This is no longer true in the presence of parity kicks: (b) shows that, after t = lOí^ec, the 
state is almost indistinguishable from the initial one and that the quantum wiggles of the Wigner function are still 
well visible. 



V. CONCLUDING REMARKS ON THE POSSIBLE EXPERIMENTAL APPLICATIONS 



The numerical results presented above provides a clear indication that perturbing the dynamics of an oscillator with 
an appropriate periòdic pulsing can be a highly efficient method for controlling decoherence. However we have seen that 
a significant decoherence suppression is obtaincd only when the pulsing frequeney 1/(T + To) becomes comparable 
to the cutoff frequeney of the reservoir. This fact poses some limitations on the experimental applicability of the 
proposed method. For example the parity kick method is certainly unfeasible, at least with the present technologies, 
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in the case of optical modes in cavities. In this case, realizing a parity kick is not difncult in principle since impulsivc 
phase shifts can be obtained using electro-optical modulators or a dispersive interaction between the optical mode 
and a fast crossing atom jï3]. The problem here is that it is practically impossible to make these kicks sufhciently 
freqüent. In fact, we can assume, optimistically, that the frequency cutoff uj Ci even though larger, is of the order of 
the cavity mode frequency ljq (this is reasonable, since the relevant bath modes are those nearly resonant with the 
frequency u)q). This means lü c /2tt ~ 10 14 - 10 15 Hz and it is evident such high vàlues for \/(T + r ) are unrealistic 
both for electro-optical modulators and for dispersive atom-cavity mode interactions. A similar situation holds for 
an electromagnètic mode in a microwave cavity as that studied in the Schròdinger cat experiment of Brune et al. 
||. In this experiment the dispersive interaction between a Rydberg atom and the microwave mode is used for the 
generation and the detection of the cat state and therefore the realization of 7r-phase shifts is easy. The problem again 
is to have sufficiently fast and freqüent 7r-phase shifts: in this case the mode frequency is lv = 51 Ghz and therefore 
one would have a good decohcrence suppression for T + tq ~ 1CP 10 sec. This is practically impossible because it 
implics quasi-relativistic velocities for the crossing Rydberg atom and an unrealistic very high dispersive interaction 
in order to have a 7r-phase shift. 

The situation is instead different in the case of ions trapped in harmònic traps. In this case in fact, the free 
oscillation frequency wo, and therefore the cutoff frequency lo c too, are usually much smaller and it becomes feasible 
to realize fast and freqüent 7r-phase shifts. Let us consider for example the case in which the oscillator mode of the 
preceding sections is the center-of-mass motion of a collection of trapped ions in a Paul trap, as that considered in 
the experiments at NIST in Boulder |Ï4|| . In this case the free oscillation frequency is of the order of lj /2tt ~ 10 Mhz 
and a freqüent sequence of parity kicks could be obtained in principle by appropriately pulsing at about 10 Mhz the 
static potential applied to the end segments of the rods confining the ions along the z-axis O). The duration t and 
the intensity E a of the pulses have to be tailored so to have the parity kick condition E t = nh (see Eq. (p0|)) and 
this implies having a very good control of the pulse area. 

Controlling the decoherence of the center-of-mass mode is crucial for the realization of quantum imformation 
processing with trapped ions. In fact, even though only the higher frequency vibrational modes will be used for 
quantum gate transitions involving the ions motion (as it has been already done in [pof where the deterministic 
cntanglement of the internal states of two ions has been achieved by manipulating the stretching mode), suppressing 
decoherence of the center-of-mass motion is still important because the heating of the center-of-mass motion partially 
couples also with the other vibrational stretching modes ]Ï4| ], 

The parity kicks method could be used even in the case in which the center-of-mass mode itself is used as quantum 
bus for the realization of quantum gates. However in this case, the use of parity kicks cannot be applied to protect 
against decoherence all kinds of quantum gates. In fact the parity kicks tends to average to zero any term in the 
Hamiltonian which is an odd function of the bosonic operators a and , and this means that the parity kicks could 
average to zero just the system dynamics we want to protect from decoherence, as for example, some gate operations 
in an ion-trap quantum computer. For example, it is easy to see that all the gate operations involving first red or 
blue sideband transitions, which implies having linear terms in a and a' in the system Hamiltonian, as for example 
the Cirac-Zoller controlled-NOT (C-NOT) gate jÏÏ) and the C-NOT gate experimentally realized at NIST by Monroe 
et al. p^ |, tend to be averaged to zero by freqüent parity kicks. To overcome this problem, it is however sufficient 
to restrict to quantum gates based on carrier transitions, which involve functions of a^a only, and are not therefore 
affected by parity kicks, as for example the one-pulse gate proposed by Monroe et al p3[ . 

Therefore, the parity kick method presented here could be very useful to achieve a significant decoherence control in 
ion traps designed for quantum information processing. However we have to notice that our considerations are based 
on the model of section III which is the one commonly used to describe dissipation, even though the actual sources of 
decoherence in ion traps have not been completely identified yet [Q. For example Refs. |24|,^5|] have considered the 
effect of fluctuations in the various experimental parameters, yielding no appreciable energy dissipation in the center 
of mass motion but only off-diagonal dephasing. 
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FIG. 1. Time evolution of the decoherence function rj{t) (see Eqs. (|3^) and ( 
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FIG. 2. Decoherence function rj(t) at a fixed time t versus u c T /2n: a) refers to í = 0.5/7 and b) refers to í = I/7. 
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FIG. 3. (a) Wigner function of the initial odd cat state, \tp) = N-(\ao) — \ — ao}), «o = v5; (b) Wigner function of the same 
cat state evolved for a time t = I/7 (t = lOí^ec), in the presence of parity kicks (ui c T = 3.125); (c) Wigner function of the same 
state after a time t — I/7, but evolved in absence of kicks. 
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